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Abstract 

In the present paper we discuss various aspects of the 't Hooft model for two- 
dimensional QCD in the limit of infinite number of colours in the Coulomb gauge. 
The properties of mesonic excitations are addressed, with special attention paid to 
the pionic one. The twofold role of the pion is discussed: being a genuine qq state it 
is also a Goldstone boson of two-dimensional QCD. In particular, it is demonstrated 
explicitly how the soft-pion theorems are satisfied. It is pointed out that the Coulomb- 
gauge choice seems to be indispensable in studies of hadronic observables with the 
pions involved. 
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1 Introduction 

Two-dimensional QCD in the limit Nc — > oo (the 't Hooft model was first considered in 
70-th. Since then the 't Hooft model is widely used as a toy laboratory for studies of various 
aspects of strong interactions. The theory is relatively simple, as there are no transverse 
gluons in two dimensions; moreover, in the large-Ac- limit only planar graphs are to be 
summed up, and the theory is exactly solvable. Nevertheless, this is a truly relativistic field 
theory which does have a nontrivial content, resembling in such a way realistic QCD4 case. 
Indeed, 

• the theory exhibits confinement and it is possible to demonstrate explicitly the exis- 
tence of the discrete spectrum of the quark-antiquark bound states; 

• the Poincare invariance is maintained; 

• the chiral symmetry is spontaneously broken; 

• the Goldstone boson responsible for the chiral symmetry breaking is the qq ground 
state. 

The first point is an almost trivial statement, since the Coulomb force is confining in two 
dimensions. 

The second item is of a paramount significance for hadronic spectroscopy. It was demon- 
strated explicitly in that, if a nonabelian theory is quantized in the explicitly noncovariant 
gauge, a special care should be taken of the Lorentz invariance. The quantum Poincare al- 
gebra is closed only in the colour- singlet sector, which means that the spectrum can be 
evaluated in an arbitrary frame including, for example, the centre-of-mass frame as well as 
the infinite-momentum one. 
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The chiral issue was historically a bit controversial. The initial studies in QCD2 were 
performed in the light-cone gauge. The pioneering paper was followed by detailed studies 
of spectra and wave functions of mesons as well as hadronic interactions |3|, [|]. A bit later 
an alternative approach was suggested, based on the Coulomb gauge A\ = 0. The main 
advantage of the light-cone quantization is considerable simplification of the spectra calcu- 
lations, but straightforward analysis gives the perturbative vacuum. The more technically 
involved version || yields a nontrivial vacuum, and it appears that a nonzero quark con- 
densate exists for the massless quarks |J. The latter feature is confirmed by the sum rule 
calculations in the light-cone gauge HQ. At the hadronic level the apparent discrepancies 
are connected with a peculiar form of the pionic wave function near the chiral limit in the 
light-cone quantization, as it is discussed in 0, and will be explicitly demonstrated below. 
That is why one is forced to employ the sum rules and the Operator Product Expansion 
(OPE) to arrive at reliable results in the pionic physics and the vacuum structure ||. On the 
contrary, the choice of the Coulomb gauge does not lead to drastic singularities and enables 
treating the pions on the same footing as other mesons. The conceptual difficulties of the 
light-cone quantization were resolved in the formulation on finite intervals |§, where the 
equivalence of both versions was demonstrated explicitly, clarifying the relationship between 
the light-cone and the equal-time quantization. 

Since then a lot of work in two-dimensional QCD was done, employing the light-cone 
gauge. Among the questions discussed are the general properties of the OPE [|, || and 
heavy quark expansion and duality [011. The calculations of spectra were performed in the 



framework of the so-called discretized light-cone quantization beyond the Nc — > 00 limit 
n|. A separate fascinating issue is the studies of QCD 2 with adjoint fermions flT2|. In the 



present paper we discuss the properties of vacuum and low-lying mesonic states built of light 
quarks in the Coulomb gauge, with special attention paid to the chiral issues of the theory. 

Quantization on the light-cone allows one to establish an obvious connection with the dy- 
namics of the parton model, while quantization on the ordinary time hypersurface is natural 
for another important branch of phenomenology, the constituent quark model. Indeed, the 
spectrum of QCD 2 is discrete, and Nq — > 00 limit suppresses additional quark pair creation, 
so that the 't Hooft model is nothing but a constituent quark model exactly derived from a 
nontrivial relativistic quantum field theory. 

In the constituent quark models the confinement is usually modeled by a potential force. 
Then the gross features of the light quarkonia spectra and decays are described surprisingly 
well with an exception of the pion. Since one cannot include chiral symmetry breaking (CSB) 
effects into the constituent picture, there is no hope to reproduce the pion as the Goldstone 
boson, and soft pion theorems cannot be satisfied within any naive quark model picture. 

The CSB phenomena follow from the most general symmetry considerations and have 
nothing to do with the particular mechanism of the confinement. One possible way to include 
the soft pions into the quark model is to introduce quarks and pions on equal footing, as it 
is done in the chiral quark models (see, e.g., |13|] and references therein). Nevertheless, there 
are some results from the lattice simulations demonstrating that the confinement and 
the CSB phenomena are present in the confining phase and disappear in unison above the 
deconfinement temperature; so both phenomena are interrelated dynamically. 

A model was suggested many years ago JT5[, which connects the confinement and CSB 



(see also fll6| , where similar ideas were employed). The main ingredient of this model is 
an instantaneous three-dimensional oscillator confining force. Such a force does not follow 
from QCD and, in addition, there is no gauge and Lorentz invariance in this approach. An 
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important development was suggested in |L7|], where the QCD vacuum is parametrized by a 
set of gauge and Lorentz invariant nonperturbative gluonic correlators which are responsible 
for both, the area law and the chiral condensate formation. The quark model which follows 
from such an approach should be able to reproduce, inter alia, all pion properties. In this 
regard it is instructive to study an exactly solvable theory with confinement and CSB, and 
the 't Hooft model is a perfect toy laboratory for such studies. 

Before proceeding further we would like to note that the large- Nq limit is essential in 
establishing the chiral properties of QCD 2 ||18|| . The Coleman theorem |19| prohibits CSB 
for any finite number of degrees of freedom in a two-dimensional theory. It means that limits 
Nc — > oo and m q — > are not interchangeable. There is no contradiction with the Coleman 
theorem if one considers the weak coupling regime where m q ^> g ~ l/y/Nc, i-e., the limit 
of infinite number of colours is taken first (see |J for the detailed discussion of this issue 
as well as of the other phase of the theory which corresponds to the strong coupling regime 
m q < g). 

The theory is given by the Lagrangian 

L(x) = -If^(x)F^{x) + q(x)(W - m)q(x), (1.1) 

D = (d,- igAffl-rn, 
and the large- Nc limit implies that 

9 2 N C , n x 

7 = — — -> const. (1.2) 

47T N c -^oo 

Our convention for 7 matrices is 

70 = (3 = a 3 , 71 = ia 2 , 75 = a = 7071 = o x . (1.3) 
The theory is gauged by the condition 

A 1 (x o ,x)=0, (1.4) 

so that the only nontrivial correlator of the gluonic fields, the gluon propagator, takes the 
form 

Doi(x -y ,x-y) = Df x = (x -y 0i x-y) = 0, 

^00(^0 -Vo,x-y) = -^ ah \x ~ y\S(x - y ), (1.5) 

and the infrared singularity is regularized by the principal value prescription, i.e., by means 
of an appropriate number of subtractions, for example, 

= J j^-^(F(k)-F(p)- F\p)(k-p) -...), (1.6) 

yielding a linear confinement for the interquark interaction mediated by the two-dimensional 
gluon. 
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(a) (b) 




Figure 1: The planar (figure (a)) and the nonplanar, suppressed by Nc, (figure (b)) diagrams. 



The paper is organized as follows. In Section 2 a Hamiltonian approach to QCD 2 in 
the Coulomb gauge is developed. The bosonization of the model is performed explicitly 
in the large- Nc limit, and a generalized Bogoliubov-Valatin transformation is employed to 
construct the composite operators creating and annihilating bosons. The pion wave function 
is found explicitly and the chiral properties of the theory are discussed. In Section 3 a matrix 
formalism is presented, a matrix bound-state equation is derived with special attention paid 
to the problem of the truncated Hilbert space and Hermiticity. We clarify the role of the 
backward motion of the qq pair in the meson. The current conservation and the Ward 
identities are discussed as well as the pionic vertex. Section 4 is devoted to strong hadronic 
decays, where we check the low-energy theorems, including the "Adler zero" selfconsistency 
condition. Concluding remarks are given in Section 5. 



2 Hamiltonian approach 

As is known from the pioneering work |I| , the physical degrees of freedom in two-dimensional 
QCD in the weak-coupling regime are qq mesons. It would be quite natural then to reformu- 
late the model entirely in terms of compound mesonic states, introducing a nonperturbative 
vacuum and creation and annihilation operators for the mesons, as excitations over this 
vacuum. The most natural framework for such a task is the Hamiltonian approach to the 
model, which we develop in this section.^] Note that it is possible in a selfconsistent form 
due to an instantaneous type of the interaction induced by (|1.5|) . 

This section is organized as follows. In subsection |2TT| we start from the Hamiltonian of 
the model in the Coulomb gauge, introduce dressed quark fields and, following 0, derive a 
gap equation (also known as mass-gap equation). We discuss the numerical solution to the 
gap equation found in ||. Subsection [27| is devoted to investigation of the vacuum energy 
and identification of the chirally-symmetric and nonsymmetric vacua. In subsection [2.3| we 



introduce two-particle operators and perform a generalized Bololiubov transformation to 



diagonalize the Hamiltonian of the model in the mesonic sector. In subsection ^]4| we discuss 



properties and solutions of the bound-state equation, which appears as a consequence of 



lr The Hamiltonian approach to two-dimensional QCD in the light-cone gauge was developed in 
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the second Bogoliubov transformation. A special solution to the bound-state equation, 
the massless chiral pion, is found analytically and investigated in subsection |2l|. An issue 
connected to the locality and the Lorentz nature of confining interaction in two-dimensional 
QCD is discussed in subsections [2.6| and [2.7| . In conclusion we turn to the chiral properties 
of the model and this is the subject of subsection [2TB. 



2.1 Dressed quarks and the mass-gap equation 



Starting from the Lagrangian and following the standard rules one arrives at the Hamil- 



tonian of the theory in the form 



H = J dxq + (t, x) (~*75T^: + ™7o^) q(t, x) 



J dxdyq + (t,x)t a q(t,x)q + (t,y)t a q(t,y) 1 -^^. (2.1) 



Note that only four-quark interaction enters the Hamiltonian ( |2.1| ) and this is a reflection 
of the fact stated in the Introduction that the only nontrivial gluonic correlator is the gluonic 
propagator (|1.5| ). In four- dimensional QCD, possessing a much more complicated vacuum 
structure, correlators of all orders should appear, which gives rise to an infinite number of 
terms with four-quark, six-quark, etc interactions in the Hamiltonian. In the meantime, it 
is reasonable to truncate the QCD 4 Hamiltonian at the four-quark interaction level, which 
corresponds to the bilocal approximation, when only (A A) correlator is left.0 The interested 
reader can find a detailed review of the given approach in ||22|| . 

The "dressed" quark field q(x) in ( |2.1|) is defined as follows f2| 

q at (t,x) = J ^ e lkx lb a (k,t) Ul (k) + d+(-k,tM-k)}, (2.2) 

& a (f,AO|0) = d a (t,-AO|0> = 0, b+(t,k)\0) = \q), d+(t,-k)\0) = \q), (2.3) 
{b a (t,p)b+(t, q)} = {d a (t, -p)d+(t, -q)} = 2n5(p - q)5 af3 , (2.4) 

u(k) = T(k) Q , v(-k) = T{k) Q , (2.5) 

T{k) = e -^ {k)l \ 

where 9 plays the role of the Bogoliubov- Valatin angle. Greek and Latin letters denote 
colour and spinor indices, respectively. Strong interaction between quarks implies that the 
true vacuum state is described by a nontrivial 8, whereas excitations over it bring positive 
contribution to the energy. 

In what follows we shall omit the explicit dependence of operators on time. It can be 
easily restored at any intermediate step, thus giving, for example, 

b a (t,p) = b a (p)e-^\ d a {t, -p) = d a {-p)e iE ^, (2.6) 

where E(jp) is the dispersive law of the dressed quark (to be defined later). 

2 These are irreducible averages (cumulants) meant here, which are defined as ((O)) = (O), ((O1O2)) = 
(O1O2) — (O1KO2), etc Due to the Lorentz and the colour invariances of the QCD vacuum (A®) = 
and, hence, {{A a A b v )) = (A a A h v ). 
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The Hamiltonian ( |2.1| ) normally ordered in the basis ( j2.3| ) splits into three parts (L is 
the one-dimensional volume), 



H = LN C £ V + :H 2 : + :H± 



where 



£ v = J T^Tr | (75J) + m7o) A-(p) + 



7 



dk 



An J (p— k) 2 



A+(A;)A_(p) 



(2.7) 
(2.8) 



is the vacuum energy density, 



H 2 := J dx: qf{x) f-^TsJ^ + ™7o^) ftfc) : - 



y dxdy- 



x-y\ 



dk:qf{x) [A + (k)-A4k)} qi (y):e ik ^ 



is quadratic in quark fields, whereas the : H± : part contains four of them, 

2 1 I 

: #4 := ~- J dxdy : q + {x)t a q{x)q + {y)t a q{y) : H^l. 



(2.9) 



(2.10) 



We have introduced projectors onto positive and negative states for convenience: 

A±(p) = T(p)^^T + (p). 



One comment on the role played by different parts of the Hamiltonian 



(2.11) 
is in order. 

The first term in (|2.7j ) defines the energy of the vacuum which is to be minimized. The 
: H2 : part describes "dressing" of quarks, so that an alternative approach to minimizing S v 
is the requirement that : H 2 : be diagonal in terms of the quark creation and annihilation 
operators, or, equivalently, that the anomalous Bogoliubov term be absent. No matter which 
way is used, the resulting equations for the Bogoliubov- Valatin angle 9 and the dispersive 
law for the dressed quarks read [0: 



E(p) cos 9{p) = m + lj + — cos 9(k) 

J [p — k) 

E(p) sin 9(p) =p + 2/7—^)2 sin0 ( k )i 
which can be reformulated in the form of the gap equation for the angle 9, 



dk 



(2.12) 



p cos 9(p) — m sin 9{p) = — j- 



dk 



sm[9(p) - 9(k)). 



(p — k) 2 

If a solution for 9 is known, E(p) can be easily found from the relation 
E(p) = mcos8(p) +p sin 8(p) + — 4 — -qos[9{p) — 0{k)\. 



(2.13) 



(2.14) 



The gap equation ( |2.13|) is subject to numerical investigation which was performed in ||. 
The results for 9(p) and E(p) are given in Fig.|2|. 

Notice several important properties of the functions 9 and E: 
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Figure 2: Numerical solutions for the Bogoliubov angle 9(p) and the dressed quark dispersive 
law E(p) for different values of the quark mass. The plots are taken from ||, where x comes 
from the change of variable p = tan(x), all dimensional quantities are given in the proper 
units of (27) 1 / 2 . 



9(p) is odd, whereas E(p) is even, 

9(-p) = -9(p), E(-p) = E(p); (2.15) 

solution for 9(p) remains nontrivial even in the chiral limit, m = 0, and so does the 
solution for the dispersive law (see FigfJ); 

the function E(p) is not positively defined, as it might be naively expected (see FigfJ). 



The formal reason for the last property comes from the regularization prescription (|L 
which defines the way of subtracting an infinite self-energy constant. Thus a divergent 
integral of positively defined functions might lead to a negative result after the proper reg- 
ularization. 

One can easily check that 

0(p^O)~p, 0'(O)>O, (2.16) 

and 

E(0) =m + lj d ^ C os9(k) = m -lJ^ sin 2 B M « m - ^V(O), (2.17) 

so that E(0) becomes negative for m smaller than some critical value m c . This does not 
lead to a disaster as, according to findings of @, only the colour-singlet sector of the theory 
is Lorentz covariant, whereas the quark sector is not colour-singlet and, hence, it may not 
be Lorentz covariant. Nevertheless, one has to be especially careful shifting the poles in 
the quark propagator, paying attention to the sign of the product of E(p) and the infinitely 
small positive constant e 0. 
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Note that a simple analytical solution to the gap equation ( |2.13| ) in the chiral limit was 
found in 0, which reads[] 

9(p) = ivrsign(p), E(p) = \p\ - P^- (2.18) 
2 \p\ 



where the symbol P stands for the principal value prescription. It is clearly seen from ( [2.18D 
that for this solution E(p) is not positively defined either. 



The two solutions to the gap equation in the chiral limit, the one given by (|2.18|) and 
the other, depicted in Fig.^|, define two different phases of the theory. The chiral symmetry 
is broken in the latter case only, whereas in the former one it remains unbroken. 



It was demonstrated in that the analytical solution ( 2.18| ) possesses infinite energy 



compared to the numerical chirally-nonsymmetric one given in Fig.|2|. Thus the chiral sym- 
metry is never restored in the 't Hooft model and there is only one phase of the theory with 
chiral symmetry spontaneously broken (see the next subsection for the details). Such way, 



the solution ( 2.18 ) is unphysical, but we still prefer to keep it to exemplify some statements 
concerning the chiral properties of the model. In the next subsection we discuss this issue 
in more detail. 

Performing all necessary calculations with the precautions discussed above, one arrives 
at the Hamiltonian of the model which is diagonal with respect to the dressed quarks basis. 
The contribution of the : : part of the Hamiltonian ([2.7D is suppressed by an extra factor 
1 / V^c and thus it can be neglected in the single-quark sector of the theory; 

: H := LN C £ V + J ^-E(k) {b + a {k)b a {k) + d+{-k)d a (-h)} . (2.19) 

Note that the result ( |2.19| ) itself has practically no value as it deals with the gauge- 
and Lorentz-noncovariant sector. The most important result of this subsection is the gap 
equation ( p. 13 ) and the numerical solution to it depicted in Fig.[|. They will be intensively 



used in what follows. 



2.2 The vacuum energy. Chirally-symmetric and nonsymmetric 
vacua 



Let us return to the gap equation (|2.13j ) and discuss an alternative way of its derivation 



namely, the minimization of the vacuum energy ( |2.8| ) [p3fl . We consider the case of massless 
quarks, m = 0. 

It is convenient to introduce an excess of the vacuum energy density for the theory with 
interaction over the free-theory one, 

A£ v [0\=£ v [O]-£ v [O fM ] = -J ^(psin^(p) - - 2- j cos[0(p) - 0{k)\, (2.20) 

where 9f ree (p) = |sign(p) corresponds to the free (7 = 0) massless theory. 
The gap equation, following from the minimization procedure, 

r 

A£ v [9} = 0, (2.21) 



59(p) 



3 T 



In fact, any odd function 9{p) with an arbitrary number of jumps from tt/2 to ~ir/2 and back also 



satisfies the gap equation (2.13) M. Meanwhile solutions of such a type do not reduce to the free theory 



when 7 tends to zero and they will not be discussed. 
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readily reproduces equation ( |2.13| ). 

To ensure that the solution to equation ( |2.13|) indeed minimizes the vacuum energy we 
use the following approach Let 0{p) be the solution to fl2.21p corresponding to the 

minimum of AS V . Then 8(p/A) stretched with an arbitrary parameter A should enlarge the 
energy (|2.20|). Naive dimensional analysis demonstrates that AS V scales with A as 



-dA 2 - 7 C 2 , 



(2.22) 



with C\ t 2 being positive constants. Then the stable solution is provided by Aq minimizing 
the energy ( |2.22j ), i.e., by Aq = 0, which corresponds either to the free massless theory, giving 
A£ v = 0, or to the analytic solution ( |2.18| ). Both solutions correspond to the preserved chiral 
symmetry. Thus one arrives at a discouraging conclusion that there is no nontrivial chirally- 
nonsymmetric solution to the gap equation fl2.13|) . In the meantime, the naive analysis 
performed above completely ignores the fact that the vacuum energy ( |2.20| ) is logarithmically 
infrared divergent due to the second term. Introducing a cut-off and repeating the same steps, 
one can conclude that the correct dependence of the vacuum energy on the scale parameter 
A is as follows: 

A£ v = \dA 2 - 7 C 2 In A + 7 C 3 (2.23) 



instead of the naive form ( |2.22| ). The constant C3 is proportional to the logarithm of the 
cut-off and can be removed by an infinite renormalization. 

Note that the second term in ( |2.20| ) cannot be made convergent both, in the infrared and 



in the ultraviolet simultaneously. Indeed, one can remove the infrared divergence in ( 2.2C ), 
subtracting the energy corresponding to the solution ( |2.18| ) instead of the free one. Then 
the resulting integral appears ultraviolet logarithmically divergent and leads to the same 
relation (|2.23|) but with the infinite constant C3 containing the logarithm of the ultraviolet 
cut-off. 

The function ( 2.23| ) always has a minimum at 



A 




[2.24) 



which corresponds to the nontrivial solution of the gap equation found numerically in || 
and depicted in Fig|J. In the meantime, from the form of the function ( |2.23|) one can see the 
logarithmic growth of the energy in approaching the solution ( 2.180 , which, as was discussed 
above, corresponds to A = 0. 

Thus ones arrives at the conclusion, already mentioned above, that the vacuum energy 
corresponding to the chirally-symmetric solution Q2.18p is infinite compared to that for the 
chirally-nonsymmetric one depicted in Fig.0] and, hence, no phase transition of the chiral 
symmetry restoration is possible in the 't Hooft model. 

It is instructive to note that the very fact of the existence of the chirally-nonsymmetric 
nontrivial solution to the gap equation ( 2.13Q is yet another consequence of the infrared 
behaviour of the 't Hooft model discussed above. Indeed, the gap equation (|2.13|) was derived 



4 Once the logarithmically divergent term in ( 2.23| ) is proportional to the coupling constant 7, then there 
is no problem with the free limit of the theory, which also corresponds to Aq — 0. Indeed, when 7 tends to 
zero the logarithmic term in (2.23) disappears, so that, as defined by ( [2.20| ), the vacuum energy of the free 
theory is zero. 
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in neglection of all effects of the fermionic vacuum polarization (creation and annihilation 
fermionic operators introduced in ( |2.3|) correspond to the so-called BCS approximation). In 
the meantime, the chiral symmetry can be spontaneously broken in other two-dimensional 
models for QCD, like the Gross- Neveu one p3, p5|, but in order to have a nonzero chiral 



condensate one has to go beyond BCS level, summing up fermionic bubbles, whereas the 
equation similar to ( |2.13| ) has only a trivial chir ally- symmetric solution, giving (qq) = 0. 



2.3 Generalized Bogoliubov transformation and mesonic compo- 
und states 



In the previous subsection the first two terms of the Hamiltonian ( |2.7Tj were considered. 
Let us turn to the third part, : if 4 :. With substituted solution for the dressed quarks, 
it describes interaction between them, which leads to formation of the qq bound states - 
mesons. In this subsection we are to diagonalize the Hamiltonian ( |2.7| ) in the colour-singlet 
mesonic sector of the theory. To this end we introduce compound operators which act on 
colourless pairs of quarks and antiquarks [0, p6fl: 



B(p,p>) = lb+(p)b a (p'), D(p,p') = ld+(-p)d a (-p f ), 



'N, 



c 



'N 



c 



M (p,p>) = ld a (-p)b a (p>), M+(p,p>) = 1- b+(p')d+(-p), 
\JN C ^N c 

with the commutation relations being 

[M(p,p')M+(q,q')} = (2n) 2 6(p' - q')5{p - q)- 
2ir 

- -= {D(q,p)6(p' - q') + B(q',p')S(p - q)} -> 



[B(p,p f )B(q,q')) 



2tt 



(2ir) 2 5(p'-q')6(p-q), 

(B(p, q')5{p' -q)- B(q,p')6{p - q')) -> 

I\ C — >00 



[D(p,p')D(q, q')} = (D(p, q')8{p' - q) - D(q,p')S(p - q')) 



>N C 

[B(p,p')M(q,q r )}-- 
[B(p,p')M+(q,q')} 
[D(p,p')M(q,q')]-- 
[D(p,p')M + (q,q')} 



N c ^oc 



2tt 



'N c 
2?r 



■M(q,p')S(p-q') 



'N c 
2tt 



M+(q,p)5(p' - q') 



'N c 
2n 



M(p', q')5{p - q) 



TVp— >oo 



N c ^oo 



M + (p, q')5(p' - q) 



(2.25) 



(2.26) 



iV c — >oo 
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With the new operators substituted, the Hamiltonian ( |2.7| ) takes the form 

dk 



I fun 

H = LN C S V + ^JN C J —E(k){B(k, k) + D(k, k)} 
dp dk dQ 



(27r) 2 (p - ky 



6(p) - 6(k) . 6(Q-p)-6(Q-k) 
2 cos sin 



x 



(2.27) 



2 2 
{M + (p,p - Q)D{k -Q,k) + M + {p - Q,p)B{k - Q, k) 

-B{p,p- Q)M(k -Q,k)- D{p,p- Q)M(k, k - Q)} 

6(p) - 9(k) e(Q-p)-9(Q-k) 

+ cos cos 

2 2 

x {B{p - Q,p)B{k, k-Q) + D{p,p- Q)D{k - Q, k) 

+ M + (p - Q,p)M{k -Q,k) + M + (p,p- Q)M(k, k - Q)} 

. 6(p) - 6(k) . 6(Q-p)-6(Q-k) 

+ sin sin 

2 2 

x {B(p,p- Q)D(k, k-Q) + B{p- Q,p)D(k - Q, k) 

+M(p,p- Q)M(k -Q,k) + M + (p - Q,p)M + (k, k - Q)} 

where : H2 : and : H4 : terms should have the same order in powers of Nc and thus act 
together as opposed to the one-body sector, where : H4 : was suppressed as 1/ \fN c . 

A crucial step we are to perform now is to note, that in the mesonic sector of the theory 
one cannot create and annihilate isolated quarks and antiquarks. Only colourless qq pairs 
can appear, so that, creating a quark, we have to create an accompanying antiquark and, 
vice versa, for each created antiquark we have an extra quark. Thus the operators ( |2.25| ) 
cannot be independent. Indeed, it is easy to check that the substitution 



B(p,p>) = ~hl -^M + (q",p)M(q",p'), 
J N„ J Ztt 

dq" 



D(p,p') 



'N c 
1 



(2.28) 



2tt 



M + (p, q")M{p',q'' 



satisfies the commutation relations ( |2.26| ). Now one can neglect a number of terms in the 
Hamiltonian (|2.27|) and to arrive at a simplified expression, 



H = LN C S V + 

_ 7 
2 



dQdp 
dk 



{E{p) + E(Q - p))M + {p - Q,p)M(p - Q,p) 



(p - ky 



{2C(p, k, Q)M + (p - Q,p)M(k - Q, k) 



(2.29) 



+S(p, k, Q) (M(p,p- Q)M(k -Q,k) + M + {p,p- Q)M + {k - Q, k))} 
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where 

r< i,n^ 9 ^ ~ 9 ^ 9(Q-p)-9(Q-k) 
C{p, k, Q) = cos cos , 

(2.30) 

q , . n , - HP) ~ m . 9(Q-p)-9(Q-k) 
S{p, k, Q) = sin sin . 

The operators M + and M create and annihilate quark-antiquark pairs, which are not 
mesons yet since Hamiltonian fl2.29| ) is not diagonal in terms of these operators. 
Symbolically the operator structure of the Hamiltonian ( p. 29 ), 

H ~ H + AM+M + -B(M + M + + MM), (2.31) 

resembles the one appearing in the theory of Bose gas with interaction, where the last term 
on the r.h.s. describes interaction between particles and leads to appearing of quasiparticles 
diagonalizing the Hamiltonian. Thus fl2.29|) is subject to another Bogoliubov transforma- 
tion. Once operators M + and M obey the Bose statistics, then the general form of the 
transformation is expected to be 

j m = uM + vM+ 

\ m+ = uM+ + vM, 1 6 > 

with u and v obeying the standard bosonic condition, 

u 2 -v 2 = l. (2.33) 
Of course, one needs to generalize transformation (|2.32| ), ( |2.33|) in order to take into 



account the nonlocality of the involved objects. Such a generalization takes the form [^6 

dq 



m n (Q) = I p- {M(q - Q, q) V n + (q, Q) + M + {q, q - Q)<p»(q, Q)} 



(2.34) 



2tt 



M(p-P,p)= E {m n (P)^l( P ,P)-m+(P)<f n _i P ,P)} 

n>0 L ' 



(2.35) 



M+(p-P,p)= E {m+(P) V l(p,P)-m n (P) ¥ > n _(p,P)\, 

where (p± stand for Bogoliubov-like functions u and v, so it is not surprise that they obey 
the following orthonormality and completeness conditions, which are nothing but the gener- 
alization of ( |2.33 ): 



nmi 



^(vI(p,Q)^(p,Q)-V-(p,Q)v-(p,Q)) = 4 
^(^«Q)-^(p,QK(p,0)) = o, 

oo 

£ {fP+ip, Q)<P+% Q) - Q)<P-(K Q)) = 2n5(p-k) 

oo 

J2^ n + (p,Q)^(k,Q)-^{p,Q)^ n + (k,Q)) = 0. 



(2.36) 



(2.37) 



n=0 



13 



m=0.18 
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0.16 - 



0.12 
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0.04 - 



J 0.00 
1.5 0.0 




Figure 3: The ip + and ip_ components of the ground-state wave function in the meson rest 
frame for different masses of the quark taken from ||. Note that (p± are even functions of 
p for the ground state. The variable x comes from the change p = tan(x), all dimensional 
quantities are given in the proper units of (27) x l 2 . 



The functions (p± play the role of the meson wave functions, moreover, one can easily 
establish the physical meaning of both. Namely, ip 7 ] describes the motion forward in time 
of the qq pair inside meson, whereas ^p n _ is responsible for its motion backward in time. We 
shall return to this issue later on when discussing the properties of the bound-state equation. 

The physical meaning of the transformation ([2.34 ) is quite obvious. Indeed, in the theory 
with a nontrivial vacuum there are two ways to produce a quark-ant iquark bound state. 
The first way, which works no matter if the vacuum is trivial or not, is to create the qq pair 
directly, by means of the operator M + ~ b + d + . In the meantime, if the vacuum is nontrivial 
and contains the chiral condensate (qq) 7^ 0, then one can "borrow" a finite number of 
correlated quark-antiquark pairs from the vacuum and to remove redundant particles, using 
the annihilation operator M ~ db. The wave functions ip± describe the contributions of these 
two procedures. It follows immediately from such a consideration that, e.g., for the case of 
massive quarks should be small since the condensate of heavy quarks is suppressed by 
the large quark mass. 

It is easy to check that operators fl2.34|) obey the standard bosonic commutation relations, 



m n (Q)m+(Q')] = 2n5(Q - Q')5 n 



(2.38) 



m n (Q)m m (Q') 
and diagonalize the Hamiltonian ( |2.29| ), 



mt(Q)mi(Q>) 



H = LN C £' V + W f-Q° n (Q)m+(Q)m n (Q), 

n=0 J 271 



(2.39) 



if the wave functions <p± are solutions to the bound-state equation which we discuss in the 
next subsection. The vacuum energy in ( |2.39| ) contains extra contributions compared to that 
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in ( |2.29| ), besides that the vacuum structure itself is changed, so that the real vacuum of 
the theory, differs from |0) introduced in (|2~3| ) and they are related through a unitary 
transformation. 



2.4 



The bound-state equation and properties of the mesonic wave 
functions 



As in case of the first Bogoliubov transformation performed in subsection |2.1| , when the 



gap equation ( |2.13| ) appeared as a condition of the Hamiltonian diagonalization, the second, 
generalized, Bogoliubov transformation described in the previous subsection also leads to 
an equation defining the Bogoliubov-like functions (p±. This is nothing but the bound-state 
equation for the mesonic spectrum of the model 



f [E(p)+E(Q-p)-Q Q ]<p + (Q,p) = 

/A U 
— - — 2 [C(p, k, Q)(p+(Q, k) - S(p, k, Q)<p.(Q, . 

[E(p)+E(Q-p) + Q ]^(Q,p) = 



(2.40) 



Unfortunately, analytical investigation of the system ( |2.40|) is possible only in some spe- 
cific cases (we discuss one of them in the next subsection), so that it is rather subject to 
numerical studies. In Fig.|3| taken from paper [|[] we give numerical solutions for the ground 
state for different quark masses. It is clearly seen that the ip- component of the mesonic 
wave function is suppressed for large masses of the quark in full agreement with the con- 
sideration presented above. It is also suppressed for highly excited states f5|. Indeed, the 
case of higher excitations is very close to the quasiclassical regime, where the answers for 
the spectra coincide with those of the quantum mechanical problem of the two-dimensional 



string with quarks at the ends ||27|| , which, in turn, is reduced to a quark potential model. 

Let us demonstrate how the celebrated 't Hooft equation comes out from ( |2.40| ) [0 . In 
the above mentioned paper the Lorentz covariance was proved for the colour-singlet sector 
of the theory, so that the spectrum of mass of the qq bound states, following from ( |2.40| ), 
should not change when one performs boosts, even for the limiting case of the boost into the 
infinite-momentum frame, P — > oo. One can easily verify that in this case S(p, k, Q) — > 0, 
whereas C(p,k,Q) turns into a step-like function, so that the region of integration on the 
r.h.s. of ( |2.40| ) shrinks to a finite interval, < x < 1, with x = p/P being the share of the 
total momentum carried by the quark (the 't Hooft variable). The y?_ component of the 



5 It is straightforward to generalize the bound-state equation ( |2.40|) for the case of a many- flavour theory. 
One needs to make the following modifications: 



E(p) -> E h (p), E(Q -p)^ E h (Q-p), p± 



Cp, k, Q) - c hh („, k,Q) = c„ 8 -»/■(*) M W-1;'*W-" . 

Sfet.Q) -» S hh (p. k,Q) = sin *aM-«aW sin M^M^j) , 



where /i and ji stand for different quark flavours. 
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wave function dies out, and the system ( |2.40| ) reduces to a single equation for <p + , 

V {x)= lim J^- v+ (p,Q), (2.42) 

coinciding with the one derived by 't Hooft in [p]]. 

For further references we give here a couple of properties of the wave functions y>±. 

vl(p T P,TP) = V n ± (p,±P), (2.43) 

<pl(P-p,P)=r }n <pl(p,P), (2.44) 

and the parity of the state is (— l) Vn+1 . The latter property allows one to classify all mesonic 
states into two groups: odd states, pions, with r\ n = 1, and even states, a's, with r\ n = — 1. 
The odd and the even states follow one by one in the spectrum starting from the lowest 
state, which is odd and it is expected to be massless in the chiral limit. This state is nothing 
but the celebrated chiral pion. 

2.5 Pionic solution 



In this subsection we return to the bound-state equation ( |2.40|) and find an exact solution 



to it. Indeed, one can easily check that the set of the following two functions: 

, fliP , Q) = ^ ( cos m^Azm ± sin w-p)+<hp) } , (2 . 45) 

N 2 JQ) = ~Q (2.46) 

7T 

satisfies the system (|2.40|) if the quark mass is put to zero. This solution turns out to be 
massless and thus it is nothing but the chiral pion — the lowest negative-parity state in 
the spectrum (see the discussion at the end of the previous subsection). As one should 



anticipate, the norm of this state (|2.46|) vanishes in the pion rest frame, whereas in the 



infinite momentum frame this solution reads 




<P*(p,Q) n ~> 0, <p* + (p,Q) n -> \ — , 0<p<Q, (2.47) 



and, after an appropriate rescaling (see equation (|2.42|) ), gives f(x) = 1, i.e., turns into the 
pionic solution found in JE], [|. 

2.6 The one-particle limit and the nonpotential quark dynamics 

Let us consider a heavy-light system containing a static antiquark source placed at the origin 
and a light quark. The most straightforward way to derive the bound-state equation for this 
system is to take the one-body limit of the system ( |2.40| ). For the infinitely heavy flavour 



one has (see the footnote on the page [15]) 



E f2 (p) = m h -> oo, cos0 /2 = l, sin# /2 =0, p hh = 0, (2.48 
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then the coefficients in ( [2.40] ) become: 



C flh (p, k, Q) = cos 6fAp) 2 9h(k \ S hh (p, k, Q) = 0. (2.49) 

From now on, in this subsection, we shall suppress all flavour indices, having in mind 
that the angle 9{p) and the dispersive law E(p) are for the light quark. 
The system ( |2.40| ) reduces to a single equation pBfl , 



e<p(p) = E(pM P ) - 7 j j-^-y 2 cos %) - W y(fc), (2.50) 

where cp = <p + , e = Q — rrif 2 . 

Note that the interaction in ( |2.50| ), given by the integral term on the r.h.s., is essentially 
nonlocal and, moreover, there is no parameter (except the mass of the quark) which could 
allow one to treat this nonlocality perturbatively, expanding the cosine under the integral. 
Meanwhile, if the quark is also heavy, then the cosine reduces to unity and equation (|2.50| ) can 



be simplified even more, yielding the local linearly rising potential between the constituents, 

2 

I p 

(E (p) + j\x\) <p(x) = e<p(x), E (p) = Jp 2 + m 2 « m + h (2.51) 

v 2m 

To get a deeper insight into the structure of the interaction in the 't Hooft model let us 
use another approach to the heavy-light system. 

2.7 A heavy-light system in the modified Fock-Schwinger gauge 
and the Lorentz nature of confinement 



In this subsection we discuss another way to derive the one-particle limit (|2.50 ) of the bound 



state equation ( |2.40| ), based on the Schwinger-Dyson equation for QCD 2 in the modified 
Fock-Schwinger gauge (Balitsky gauge) |28j . The advantage of this approach is a possibility 
of its generalization to QCD 4 if an appropriate model for the QCD vacuum is used [[Hj . 

First of all, note that the Coulomb gauge condition, A\(x ,x) = 0, does not fix the 
gauge completely allowing purely time- dependent gauge transformations. To fix the residual 
invariance we impose an extra condition, A q (xq,x = 0) = 0,0 which obviously breaks trans- 
lational invariance, but turns out very convenient if an infinitely heavy (static) particle is 
involved. In case of the system containing a light quark and a static antiquark source the 
origin can be associated with the latter, which appears to play an extremely passive role in 
the interaction providing the white colour of the entire object. The Green's function of the 
static antiquark is independent of the gluonic field in this gauge, Sq(x, y\A) = Sq(x — y), 
and takes the form: 

S Q (x) = (l±^9{-x )e iMx ° + l -Z^6{x,)e' m ^ 5(x), (2.52) 

where 1 denotes unity in the colour space. The Green's function ( |2.52j ) contains also an 
infinite set of corrections of the form ^^^(^ — y) which die out in the limit M — > oo. 



In 3+1 this gauge is usually introduced via conditions Aq(xq, 0) = and xA a (xo, x) = p£ 
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The Green's function of the colourless qQ system can be written in the following form: 
S qQ (x,y) = ±-j D^D^DA^exp^-'-j 'd 2 xF« 2 + t J d 2 x^(id - m - A)^} x (2.53) 

x ^(x)Sq(x - y)ip(y), 

so that it turns out possible to integrate out the gluonic field arriving at a Schwinger- Dyson 
equation for the light-quark Green's function S(x,y), 

(id x - m)S(x, y) - (2tt) 2 7 / d 2 z lo S(x, z) lo b m (x, z)S(z, y) = 6 {2 \x - y). (2.54) 



Note that S(x,y) = j^S%(x,y) possesses all properties of the full qQ Green's function 
due to the passive role of the antiquark discussed above. Then both, one-particle (e.g., the 
chiral condensate), as well as two-particle (e.g., the spectrum of bound states) properties 
of the system can be extracted from the single function S(x,y). A special attention is to 
be payed to the gluonic propagator D 00 (x,y), which looks similarly to that in the Coulomb 
gauge ( |1.5| ) but contains extra terms breaking the translational invariance and encoding the 
light-quark interaction with the static antiquark, 

Doo(x - yo,x,y) = -^(\x -y\- \x\ - \y\)6(x - y ) = K(x,y)8(x - y ), (2.55) 
or in the momentum space, 

K( P ,q)=K il \ P ,q) + K {2 \ P ,q), 

K^(p,q) = \5(p-q), (2.56) 

K^( P ,q) = -- 2 5( P )-- 2 5(q), (2.57) 
qz pz 

where we have separated the local and the nonlocal parts. 

From now on two different strategies can be adopted, which finally lead to the same 
equation for the spectrum of the heavy-light system. The first approach is based on a 
diagrammatic technique with two different internal lines prescribed to the local and the 
nonlocal parts of the kernel P8| . It turns out that, in spite of the nonlinearity of the equation 
( |2.54|) , the two parts of the kernel can be considered separately due to very peculiar properties 



of rainbow diagrams |2"Bfl . Thus the local part defines the mass operator S(p), which 
can be naturally parametrized by means of the dressed quark dispersive law E( P ) and the 
chiral angle 6( P ), 

E(p) = [E(p) cos 9(p) - m] + 7i [E(p) sin 6(p) - p] , (2.58) 

with the system of coupled equations ( |2.12|) being the selfconsistency condition of such a 
parametrization. Then the nonlocal part eventually gives the bound-state equation (see 
[p8"l for the details). 

Here we choose the other strategy based on the Foldy-Wouthoysen transformation of 
equation ( |2.54|) pOfl . First, we rewrite this equation in the momentum space and use a 
spectral decomposition for the light-quark Green's function, 

a,*,,*) = arfta. - *») ( E tf>fal) ^' ( n 2) + £ ^2tM) , (2 . 59) 
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where the positive- and the negative-energy solutions ip^' have been introduced. To proceed 
further we assume that a Foldy-Wouthoysen operator, T F (p) = e~^ eF ^ 1 , diagonalizing 
equation (|2.54 ), exists and that the angle 9 F is the same for all n's. With such an assumption 
one has 

^\P) = ^(p)I>(p)fiY V { n\p) = <P n (p)T F (p)(% (2.60) 



1 



/|^»(PR(P) = U E^»(PK(9) = 2 ^(P-?). ( 2 - 61 ) 
so that the following relation holds true for the Green's function ( p. 59 ): 



/ 



( ^S(uj,q 1 ,q 2 ) = -m8{qi - g 2 )[cos6» jF (gi) - 71 sin #F(gi)], (2.62) 

Z7T 



where to = q 10 - q 2 o- 

The Schwinger-Dyson equation (|2.54j) reduces then to a Dirac-type equation in the Hamil- 
tonian form, 

(ap + (3m)Lp® n (p) - — J dqdk{(5 cos 9 F (q) + asva. 9 F (q)) K (p— q,k — q)(p G n (k) — e n (p° n . (2.63) 

The local part of the interaction in ( |2.63| ), generated by leads to dressing of the 

light quark described by the Bogoliubov-Valatin angle 9(p) and the dressed quark dispersive 
law E(p) obeying the system ( 2.12 ). Therefore, one comes to the conclusion that the Foldy 
angle 9p{p) coincides with the Bogoliubov-Valatin one, 

9 F (p) = 9(p). (2.64) 

The nonlocal part of the interaction in ( |2.63| ), which stems from K^ 2 \ is also diagonalized 
then, so that one ends with a Schrodinger-type equation, 

enVlip) = E(pM(p) - 7 / cos e{jP) -° {k \ l{k), (2.65) 



which coincides with 

Comparing bound-state equation ( p.65|) with ( p.40| ) one finds that f° plays the role of 
the ip + component of the heavy-light system wave functions, whereas ip- vanishes due to 
presence of the infinitely massive antiquark. Thus relations ( |2.61| ) follow immediately from 
( |2.36|) and (|2 . 3 7| ) with all y9_'s put to zero. 

It is instructive to note that in the Coulomb (as well as Balitsky) gauge the 't Hoof model 
is totally defined by only one nontrivial function 9(p), solution to the gap equation (|2.13j ), 
which plays a threefold role: 

• it defines the Bogoliubov-Valatin transformation from bare to dressed quarks; 

• it gives the Foldy angle, which diagonalizes the interquark interaction in the model; 

• it entirely defines all quantities in the model, including the bound-state equation. 
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Several comments concerning equation ( [2.54 ) are in order here. The first one deals with 



the generalization of (|2.54 ) to the four- dimensional case. The attentive reader may notice 



that the only two-dimensional constituent of ( |2.54j ) is the gluonic propagator D^ix, z) taken 



in the form ( 2.55| ). The equation itself survives in the case of QCD 4 , if one has an appropriate 



form of the bilocal gluonic correlator D^ u (x, z) in the given gauge, and some arguments exist, 
why higher orders correlators, which lead to many-fermion vertices higher than four, can be 



neglected (see |31j for details). 

Another interesting issue concerning equation ( |2.54| ) is that one cannot simplify the 
interaction kernel substituting 

7oS7o -> 7oSo7o, (2.66) 



as proposed in |32j. The reason for this failure is discussed in detail in and comes 
from the fact that the real parameter defining the substitution ( |2.66| ), with the consequent 
expanding of Sq in powers of the one- dimensional momentum, is the product of the quark 
mass and the gluonic correlation length. The latter parameter defines also the radius of the 
string formed between the colour constituents in the theory. A simple dimensional analysis 
demonstrates that strings are infinitely thin in 1+1, as the system has too low dimension to 
allow them to swell. One can arrive at the same conclusion inspecting the two-dimensional 
correlator (FF) = Tr(F fil/ (x)$(x,y)F pa (y)§(y,x)), where y) is the standard parallel 
transporter along an arbitrary path between the points x and y, which provides the gauge 
invariance of the entire nonlocal object ]2T] (see also the review paper for the detailed 
description of the formalism). Using the gluonic propagator ( |1.5|) , one easily finds that (FF) 
is proportional to the two-dimensional 5-function in the configuration space, 

(FF) ~ 5 <y2 \x — y), (2.67) 

i.e., it has zero correlation length T g . Thus the product vnT g is identically zero in two- 
dimensional QCD, which makes the interaction essentially nonlocal and the quark dynamics 



becomes not potential [33, 



Finally, equation ( |2.63p answers, at least in QCD 2 , the long-standing question on the 
Lorentz nature of confinement. One should distinguish between the Lorentz structure of the 
confining interaction which is of the 70 x 70 type (see ( |2.54|) ) and the effective interaction, 
which enters the Dirac-like equation ( |2.63| ), describing the bound-state problem. As clearly 
seen from ( p. 63 ), the latter contains only effective scalar (terms ~ 70) and space vector 



(terms ~ 7071 = 75) interactions. 



2.8 The chiral properties of the model in the Hamiltonian ap- 
proach 

In this subsection we discuss the chiral properties of the model which are highly nontrivial 
and in many features resemble those of four- dimensional QCD. 

Most studies of the 't Hooft model have been performed in the light-cone gauge, A_ = 
0, which leads to a perturbative vacuum and to a simpler bound-state equation. In the 
meantime, one has to employ rather sophisticated methods to discuss the chiral limit of the 
model, when all functions and distributions become extremely singular, and the whole range 
of the 't Hooft variable definition, < x < 1, is squeezed to small intervals near the boundary 
points x = and x = 1, which deliver all nontrivial content of the theory. The Hamiltonian 
approach in the Coulomb gauge developed above is free of this drawback. Indeed, to calculate 
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a matrix element of any operator between mesonic states, the only relevant ones in the weak 
regime of the model, one is to rewrite the above-mentioned operator in terms of operators m + 
and m introduced in (|2.34Q and to use the second quantization technique to evaluate directly 



the matrix element. The result is always expressed in terms of trigonometric functions of 
the angle 8(p) and integrals of them, which can be worked out analytically in some cases, or 
treated numerically. Anyway, with the numerical solution for 9 found in || any value in the 
model appears calculable. 

Let us make just one more comment concerning the vacuum of the theory. As it was 
mentioned above, the true vacuum of the model is the mesonic one annihilated by the mesonic 
operators, 

m n (P)\n) = 0, (2.68) 

for any n and any total momentum P. By a unitary transformation this state is related to 
the quark vacuum defined by relations 



= U\0). (2.69) 

Despite of the fact that the explicit form of the operator U is unknown, the difference 
between the averages calculated with the help of the mesonic and the quark vacua turns out 
to be suppressed in the large-iVc* limit, so that it plays the role of a small correction and, 
hence, lies beyond the scope of the present paper. Thus for practical calculations one is free 
to use any of the above two vacua. 



2.8.1 The chiral condensate 

A crucial test for the chiral symmetry, to see if it is spontaneously broken or not, is the 
chiral condensate (qq) . If this average does not vanish for the vacuum state, then the whole 
tower of physical states will lack the chiral symmetry respected by the Hamiltonian of the 
theory, so that this symmetry appears spontaneously broken. For the simplest case of only 
one quark flavour the chiral symmetry breaking reads U(1)l x U(1)r — > U(l)v, and it is 
very important that the U(1)a invariance is broken spontaneously in the 't Hooft model, 
in contrast to QCD 4 where this breaking is explicit due to the axial anomaly. Indeed, the 
two-dimensional anomaly is proportional to the colour trace of the coloured object F ~ 
£fiuFfj_ u , which obviously vanishes (in the meantime, the axial anomaly does exist in the 
two-dimensional QED, known as the Schwinger model, where no colour trace should be 
taken). 

Now we are in the position to evaluate the chiral condensate for the 't Hooft model. 
Following the general approach described at the beginning of this subsection, one can find 

(qq) = (Q \q a (x)q a (x)\ Q) « (0 \q a (x)q a (x) | 0) = 

r dk N + r° 

= N c / — Tr{ 7o A„(£;)} = -— dkcos6{k). (2.70) 

J 2,71 71 J 



It is instructive to arrive at the same result using another approach - - namely, the 
definition of the condensate via the light-quark Green's function. In spite of the fact that 
the single quark Green's function is a gauge variant object and, hence, it is not physical by 
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itself, when taken with the coinciding arguments and summed up over the colours, it readily 
gives the chiral condensate^ 

(qq) = -i Tr S2(x,y). (2.71) 



Once the relation ([2.71 ) is gauge invariant, then let us choose the modified Fock-Schwing- 
er gauge discussed above. The quark Green's function, S(x,y) = j^S^(x,y), is the solution 



to equation (|2.54| ) |p8fl . Using its spectral decomposition in the coordinate space, 



S(x - y ,x,y) = -i^ { n\^ { n\yy~ t£n{x °~ Vo) ^ - Vo) + 

n 

+ t^L~ ) (x)^~\y)e l£n{xo - yo) 0(y - x ), (2.72) 

n 

^\x)= r^#(p)e*", (2-73) 

where e n are the eigenenergies defined by equation (|2.65|) , substituting the Foldy-rotated 
wave functions (|2.60|) , 



(2.74) 

V { n\p) = T + (pU-\p) = V° n (p)T + (p)^}, 



and, finally, using the simplified completeness condition ( |2.61| ) for the set {<^°(p)}, one 
reproduces the result ( |2.70| ) |28| . 



It is easily seen from the definition (|2.71 ) that diagrammatically the chiral condensate can 
be represented as a closed fermion line which begins and ends at the coinciding points. Such 
way both, the motion forward in time (the positive-energy solutions) and the one backward 
in time (the negative-energy solutions), are equally important for the condensate. At first 
glance this statement contradicts the observation made above, that all </?_ components vanish 
for the heavy-light system. Solution of this problem can be found in properties of the bound- 
state equation ( |2.40| ). Indeed, there are, in fact, two sets of solutions to the system ( |2.40| ), 



with Qo > and Qo < 0, trivially connected with one another, 

<p?(p,P) = V?L(p,P), V Z n (p,P)=<pl(p,P), (2.75) 

where positive n's numerate states with Q > 0, and negative n's are prescribed to the states 
with Qo < 0. Therefore, these are ip+ n (p,P) = (pt{p,P) (n > 0) to vanish for the heavy- 
light system, whereas the two remaining wave functions (pZ n {p,P) = (p+(p,P) describe the 
propagation of the qQ system either forward or backward in time without Zitterbewegung. 
They both contribute on equal footing to the chiral condensate (|2.71|). 



From (|2.70|) one can see that the properties of the solution for 6 to the gap equation 
fl2.13p are of paramount importance for the chiral symmetry breaking. If the integral on 
the r.h.s. in ( 2.70|) vanishes, then we are in the phase of the theory with the restored chiral 



symmetry. 



7 A special care should be taken at the both stages. Indeed, the Green's function contains a discontinuity 
at x = y, so that one should approach this limit either from the side of larger, or smaller y's. On the other 
hand, if the condensate is calculated beyond the chiral limit, then the logarithmically- ultraviolet-divergent 



perturbative contribution, proportional to the quark mass, should be subtracted from ( 2.71 ) 
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It was mentioned in subsection |2.1| that there exist two different solutions to the gap 
equation ( |2.13| ) in the chiral limit. One of them, found analytically in |J, gives \6{p)\ = tt/2 
and, hence, cos 6(p) = everywhere, so that the chiral condensate ( |2.7(J| ) vanishes for this 
solution (as was discussed above, this phase has an infinite energy and, hence, never realises). 
Luckily it is not so for the numerical solution found in |5] and depicted in Fig.^. Substituting 
it into (|2.70|) and working out the integral numerically one finds: 



(qq) 



m=0 



-0.29AW27 



(2.76) 



that coincides with the results found in The chiral symmetry is spontaneously broken 
in this phase of the theory and the pion, found in the subsection [2.5|, is, indeed, the corre- 
sponding Goldstone boson. Note that the chiral condensate for the 't Hooft model is known 
analytically for any value of the quark mass 0. 

Now we can return to formula ( |2.23| ) for the vacuum energy and to rewrite it using a more 
physically transparent language. Indeed, applying the same transformation, 6(p) — > 6(p/A), 



to the chiral condensate (|2.70| ), one easily finds that E = (qq) scales linearly with A, 



(2.77) 



so that the mute parameter A can be changed for the chiral condensate and relation fl2.23|) 
can be written as 



E 
S 



; -In 



+ lC 3 



3; 



(2.78) 



where the minimum of the vacuum energy is reached for E = E given by (|2.76|) . If the quark 
mass does not vanish, then the vacuum energy density (|2.78|) acquires an extra contribution, 



AS V -> A£ v + m 



N c 



E 



(2.79) 



which explicitly breaks the invariance of AS V with respect to the change E — ► — E. Then the 
lowest (pionic) excitation over the vacuum with the wrong sign of the condensate acquires 
an imaginary mass and becomes the tachyon as it follows immediately from the Gell-Mann- 
Oakes-Renner relation (see equation (|3.20|) ). 

Minimization procedure for the vacuum energy in presence of the mass term (|2.79 ) leads 
to a more complicated equation for E, so that Eq does not provide the minimum anymore. 



2.8.2 The pion decay constant 

In this subsection we derive the decay constant f n for the chiral pion. Using the standard 
definition for it, 

(0| 4(x)\n(Q)) = UQ,^g= o , (2-80) 



i The corresponding result from || reads 



(qq)m=> 







v/6 ' 



that numerically coincides with (2.76), thus giving evidence that various integrals of 9 can be found not only 
numerically, but also in the form of irrational numbers. 
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where 

4{x) = q{x)^ 5 q{x), (2.81) 

one can calculate the matrix element explicitly using the technique described above, so that 
the result reads 

U = J^- (2.82) 



Note that the pion decay constant ( |2.82| ) is dimensionless in the 't Hooft model, which 



drastically differs in this point from QCD 4 , where this constant is dimensional and it appears 
rather small at the hadronic scale (93MeV). Thus in the four- dimensional case the pion decay 
constant defines a new scale for the effective low-energy QCD — the chiral perturbation 
theory, which cannot be developed in the 't Hooft model in view of the dimensionlessness of 

2.8.3 Partial conservation of the axial-vector current (PCAC) 

In this subsection we derive explicitly the PCAC relation for the 't Hooft model. 

Starting from the definition of the axial- vector current ( |2.81| ), we use representation ( |2.2| ) 
for the quark fields and introduce bosonic operators ( [2.25| ) after an appropriate ordering of 
the antiquark creation and annihilation operators d + and d. Leaving only the leading terms 
in the 1/Nc expansion, one arrives at 

4(x) = N C J ^v + (-k) W v(-k)+ (2.83) 



'N, 



(27Tp 

The explicit form of the quark amplitudes u and v given by ( |2.5|) together with the 
definition of the mesonic creation and annihilation operators Q2.34j ) allows one to proceed 
further and to rewrite the components of the axial-vector current in the form 



/ r dP 00 r dr> 

Jl{x) = 2^T C / N n e-^J2( m n(P)-mt(-P)) / f-g^p, P)f n (p, P), (2.84) 
J 2tt n=0 J 2tt 

i r dP 00 r dn 

J^ x ) = 2^J—N 7T e-^J2( m n(P) + m+(-P)) J £u(p,P)g n (p,P), (2.85) 



where 



f n (p,P) = ~&l(p,P)-^(p,P)), g n (p,P) = i(^(p,P) + ^(p,P)). (2.86) 
From (|2.36|) one easily finds that 

/ P)fn(p, P) = J ^Up, P)g n (p, P) = \s n7n (2.87) 

i.e., in the chiral limit the axial- vector current couples only to pions, and it can be written 

as 



4(x) = iJ^d.J (e- Wx m n (P) + e^m+(P)) = iUd^ n (x), (2.88) 



2tt ^2P q 
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where the second-quantized wave function of the pion in the coordinate space ^f n (x) is 
introduced. 

Relation ( |2.88| ) is nothing but the celebrated partial conservation of the axial-vector 
current (PCAC), whose operator form is usually formulated as a hypothesis in QCD 4 . In 
the 't Hooft model the latter can be proved explicitly and the Hamiltonian approach to the 
model turns out the most natural environment for this task. 

It is instructive to note that the form of the pionic solution can be easily guessed even 
before the : if 4 : part of the Hamiltonian is taken into account. Indeed, if the chiral sym- 
metry is spontaneously broken, then the corresponding charge does not commute with the 
Hamiltonian, 

[Q 5 H] ^ 0, Q S = J dxJ^(x), (2.89) 

so that, if the Hamiltonian is diagonalized in the quark sector, then Q 5 contains an anomalous 
term, 

Qs ~ / ^ [h+{p)d+{ ~ P) + d (~P) b (P)] cos6 (P)' ( 2 - 9 °) 

with the coefficient cos 9{p) being right the pion wave function in the rest frame (see equations 
(OD and (pi)). 



3 Matrix approach 

In spite of evident technical advantages and physical transparency of the Hamiltonian ap- 
proach developed and exploited in the previous section, it has a number of disadvantages. 
Among those are a rather tedious algebra and not straightforward connection to the dia- 
grammatic technique which is very convenient in studies of variety of hadronic processes. 
In the present section we develop a matrix approach to the 't Hooft model, which allows 
one to simplify considerably investigations of some phenomena, e.g., this technique will be 
effectively used in the next section, where the Ward identities and the strong hadronic decays 
are discussed. 

The section is organized as follows. In subsection 3.1, following [[|, we introduce the 
matrix wave function and derive the bound-state equation for it. In subsection 3.2 we study 
properties of the matrix Hamiltonian — namely, its Hermiticity and the Hilbert space of its 
definition. Chiral properties of the 't Hooft model are the subject of the next subsection 3.3. 
We establish the Gell-Mann-Oakes-Renner relation, discuss the pionic wave function beyond 
the chiral limit, and return to the calculation of the pion decay constant. In the consequent 
subsections 3.4, 3.5, and 3.6 we derive the quark-quark scattering amplitude and Ward 
identities for the vector and axial-vector currents and find the pionic vertex, respectively. 

3.1 Matrix wave functions and matrix bound state equation 

In this subsection we briefly recall the method and the results of the paper |2| based on the 
diagrammatic approach to the theory. 

At the first step we define the mass operator E as a sum of planar diagrams (see Fig.[|) 
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Figure 4: Graphical representation of the equations for the dressed quark propagator and 
the quark mass operator. 



which contribute to the dressed quark propagator S(p lx )j^\ 

1 

p — m — + ie 



S M = 2 ~ . ... . t 3 ' 1 ) 



Note that due to the instantaneous type of the interaction the integration over ko is trivial 
and the mass operator depends only on the spatial component of the momentum. Using the 
same parametrization as in ( |2.58| ) one immediately arrives at the gap equation in the form 
( p,13|) and the definition of E(p) via 9{p) coinciding with (|2.14|) . 



As the second step a homogeneous Bethe-Salpeter equation is used, which is diagram- 
matically represented in Fig.^|and defines the spectrum of the quark- ant iquark bound-states. 
The fat lines denote the dressed quark propagators (|3.1| ) whereas the meson-quark-antiquark 
vertices are described by the function T(p,P). It is also convenient to introduce a modified 
vertex f(p M ,P M ), 

f(p„,P M ) = -i 1 S(p)T(p,P)S(p-P), (3.3) 
and a matrix wave function <3>(p, Q) defined in the standard way @, 

Hp,Q) = I ir-ripM = I ^(po-Qo,p,Q). (3.4) 

J 2,71 J ATI 

The equation corresponding to the diagrams in Fig.|5| reads 

f(p M , Qfi) = ^ J f^ d *y S (P»hor(k», QfihoSiPft - QJ, (3.5) 

or, after integrating both sides of this equation over p , introduction of the wave function $ 
according to relation Q3.4j) , and performing simple algebraic transformations, one arrives at 
the bound-state equation in the matrix form, 

Qo$Go, Q) = (75P + 7om)$Go, Q) ~ $(p, Q)(l5(Q -p) + lom) 

/A V* 
j—^ {A + (fc)$(p, Q)A_(Q - k) - A+(pMk, Q)A_(Q - p) (3.6) 



9 We use the argument p M as a shorthand notation for (po,p). If not stated explicitly, then dependence 
only on the one-dimensional spatial momentum is meant, like in S(p). 



26 



Figure 5: Graphical representation for the Bethe-Salpeter equation (|3.5|). 



-A_(*)$(p, Q)A + (Q -k) + A_(p)$(*, Q)A + (Q - P )} 
where we used projectors ( [2.11 ) and the matrix wave function is parametrized as 



9(p,Q) = T(p) ( l±^ 75¥ , + ( p ,g) + 1_^ 75 ^_( P) Q) ) T + (Q-p). (3.7) 



When written in components, equation ( |3.6| ) readily reproduces the bound-state equation 
( [2.40|) for the functions ip± derived earlier in the framework of the Hamiltonian approach to 



the model. 

In conclusion we give the connection between the meson-quark-antiquark vertex T(p, P) 
and the matrix wave function P): 

r(p? p) = 1 2^ 7o fr^ 70 ' r °°' p) = 7 ° r+(p ' p)70 ' (3 - 8) 

where the vertices T and f stand for the incoming and the outgoing mesons, respectively. 
It is easy to check that, with such a definition, relation ( |3.4| ) is satisfied automatically (see 
also Appendix A for the properties of T and f ). 

3.2 Truncated Hilbert space and the problem of Hermiticity 

In this subsection we discuss properties of the matrix bound-state equation (|3.(j|), but let us 



make a comment concerning its scalar form (|2.40 ) first. As stated before, the norm of the 



wave functions ip± Q2.36|) is defined in an unusual way. Indeed, the sign minus between its 
"+" and "— " parts appears quite naturally in the Bogoliubov-like approach developed above, 
but it looks somewhat surprisingly in the context of the standard Hamiltonian technique. 
Besides, it is easy to check that, if the matrix bound-state equation ( |3.6| ) is written in the 
Schrodinger-like form, 

then the corresponding Hamiltonian IK appears non- Hermit ian. The following two questions 
should be discussed in this connection: i) what is the reason for this, and ii) whether this 
does not lead to a disaster and the eigenenergies of this equation are still real. 

The answer to the first question becomes clear if one notices that the matrix wave function 
( |3.7| ) satisfies the following conditions: 

A + (pMp, Q)A+(Q -p) = A_(p)$(p, Q)A_(Q - p) = 0, (3.10) 
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so the phase space is truncated and the Hamiltonian H acts in a subspace, that explains 
also the distorted norm ( |2.36| ). 

In the meantime, the second question concerning the spectrum persists. In order to 



answer it, let us integrate both sides of ( |3.6| ) over p, do the same for the complex conjugated 
equation, and take an appropriate linear combination. Then one arrives at the relation 

E (Qo-QD J ^ (v+(p, Q)<p+{p, Q) - <p-(p, Q)<p-(p, Q)) = 0, (3.11) 



which immediately leads to the following two conclusions: 

Qo = Qo 1 

and 



J ^ (<fl(p, Q)<p+(p, Q) - <f-ip, Q)f-(p, Q)) = S n 



^ ^P n + (p,Q)^(p,Q)-^(p,Q)^(p,Q) 



(3.12) 



(3.13) 



0. 



It was already mentioned before (subsection "Chiral condensate" above) that solutions 
of the system (|2~40D appear in pairs: for each eigenvalue Qq with the eigenfunction (</>" , </?") 
there exists another eigenvalue, — Qq, with the eigenfunction ((pi, With this symmetry, 
equation ( |3.13|) can be rewritten in the form ( [2.36|) where only positive eigenvalues enteif^]. 
Similarly in attempts to construct the Green's function for the system ( [2.40 ) the completeness 
( |2.37|) can be derived. 

Let us introduce operators C and S: 



C(S)F(p,P) = 1 



dk 



(p-k) 



-C{S){p,k,P)F{k,P) 



(3.14) 



for an arbitrary function F(p, P) with C(p, k, P) and S(p, k, P) defined in ( |2.30|) 
Then the matrix Hamiltonian 7i can be written in the form 



H 



K-C 
-S 



s 

-K + C 



l0 (K-C) + ll S, 



(3.15) 



where K = E(p)+E(P—p) is the kinetic energy. This is the term proportional to 71 in (|3.15| ), 
which makes the Hamiltonian Ti non-Hermitian. The symmetry property of the solution with 
respect to interchange of the plus and the minus components of the wave function discussed 
above follows immediately from the fact that Ti, anticommutes with 75, so that if tp = (^ + ) 

is the eigenfunction corresponding to the eigenvalue Qo, then ip' — \ t) = 75^ is a l so a 
solution with the eigenvalue — Qq, 



(3.16) 



Moreover, the eigenstate problem ( ^.9[ ) for the operator Ti. can be formulated now in the 
form of an effective Dirac-type equation, 



l0 {K-C)+ ll S-Q ^ = 



(3.17) 



10 



If not stated explicitly, we use the symbol J^ ra for summation over positive n's only. 



28 



where, as before, 

*=(£)■ (3-18) 

Mapping of the quark-antiquark bound-states problem to the properties of the fermionic- 
type equation ( |3.17| ) may be continued, which is, however, beyond the scope of the present 
paper. 



3.3 The chiral properties of the model in the matrix approach 

In this subsection we return to the chiral properties of the 't Hooft model and discuss some 
of them in the framework of the matrix formalism. 

In Appendix A we collect formulae useful for various calculations in the suggested ap- 
proach. They are entirely based on the definition of the matrix wave function ( |3.7| ) and 
properties of the matrix bound-state equation (|3.6|) . 



3.3.1 The Gell-Mann-Oakes-Renner relation and the mass of the pion 

In order to demonstrate how the matrix approach works in practice let us derive the Gell- 
Mann-Oakes-Renner relation for the 't Hooft model. We slightly relax the chiral limit in- 
troducing a small quark mass m. The matrix bound-state equation (|3.6| ) is the main object 
of investigation now. We multiply it by 7075, take trace over spinor indices and integrate 
both sides of the resulting equation over the momentum p. A number of terms containing 
projectors A± disappears and the result reads 

Qo I ^Sp[ lol5 <S>(p,Q)} -Q I ^^p[ 7l75 $(p, Q)\ = -2m I ^Sp[ l5 ^(p,Q)}. (3.19) 

J Z7T J Z7T J Z7T 

If one uses the definition of the matrix wave function and substitutes the pion solution 
( |2.45 ) into it, then relation ( |3.19| ) simplifies even more and takes the form: 

flMl = -2m(qq), (3.20) 

in which one can easily recognize the celebrated Gell-Mann-Oakes-Renner relation [34|. This 
defines the mass of the pion near the chiral limit, 

Ml = 2m I dpcos6(p). (3.21) 
Jo 



one 



With the help of the numerical solution for 9 (see Fig.0) and the footnote at page |23 
can find: 

Ml = \j— ^-JL~mSy. (3.22) 

Note that in the case of the chirally invariant vacuum, i.e., for the analytic solution 
( |2.18|) , equation (ft.20|) is trivial as its both sides vanish simultaneously. 
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Figure 6: Numerical solution for the function A and the ratio of the subleading term in ( |3.24|) 
to the leading one for m = 0.18. As before, variable x comes from the change p = tan(x), 
all dimensional quantities are given in proper units of (27) x ' 2 . 



3.3.2 The pionic solution beyond the chiral limit 



With the pion mass (|3.22j ) in hands we are in the position to go slightly beyond the chiral 
limit and to study the pionic solution in the rest frame. A simple analysis demonstrates that 
in the pion rest frame the functions g^ip, P) and f w (p, P) defined in Q2.86Q behave like 



g n (p, P = 0) ~ -±= cos 9{p) + O (Ml' 2 ) , U(p, P = 0) ~ O {^M^j . 



(3.23) 



Therefore, the leading correction to the pion wave function comes from f 1T {p 1 P) and 
<p±(p, P = 0) can be parametrized in the form (see ( [2.45D ): 



<pI(p,p = o)=n; 



-1 



cos 6{p) 

7mZ 



± JM„A(p) 



(3.24) 



where we have extracted the explicit dependence of the coefficients on the small pion mass, so 
that the unknown correction function A(p) (note that A(p) has the same parity as cos6(p), 
i.e., it is even) does not depend on M n . The dimensionless norm N n is 



Nl = Aj ^A(j9)cos%), 



(3.25) 



providing the correct normalization for tp±. 

Substituting ( |3.24|) into ( p.40|) , one arrives at a system of two equations, the first of which 
is satisfied identically due to ( |2.12| ), whereas the other one defines the correction function 
A(p), 



1 'y f dk 

- cos dip) - pA(p) sin dip) = -j , _ fc . 2 



[A(p) - A(Jfe)] cos[0(p) -0{k)\. 



(3.26) 



Equation ( |3.26| ) is subject to numerical studies which are beyond the scope of the present 
paper. 
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Figure 7: The meson- vector (figure (a)) and meson- axial- vector (figure (b)) current couplings. 



Exact numerical solutions for the pion wave functions in the rest frame for several values 
of the quark mass taken from || are given in Fig.|3], so that the function A(p) can be extracted 
directly from these data. In Fig.|6] we give the form of the function A(p) and the ratio of the 
correction defined by A(p) to the leading term in ( |3.24j ) for m = 0.18 (in units of 1/27). One 
can see from the right plot in Fig.|6] that for the given value of the quark mass the correction 
is of order one third at largest and decreases with the growth of the argument. 



3.3.3 The pion decay constant revisited 

In this subsection we give another example of calculations using the matrix approach - 
namely, we calculate the pion decay constant once again, which comes now from the fish-like 
diagram depicted in Fig.[7](b). 

The matrix element written for this diagram reads 

72 

A, = -TLNc(-I) J —^Sp [S(p)T n (p, P)S(p - P) 7/j75 ] , (3.27) 



where the factor — ry / \/ Nc comes from the vertex, whereas Nc and (-1) are due to the 
fermionic loop. Working out the integral over po and using relations ( |A.7| ) and ( |A.8| ), one 
arrives at a simple formula, 



7/^[$ n (p,P)7 M 75]. (3.28) 

Then, on substituting the explicit form of the matrix wave function $ n , putting m = 0, 
and, finally, on taking the integral over p by means of the orthogonality condition (|2.36| ), 
one finds this matrix element in the chiral limit to be 



^ = Mv F vk (3 ' 29) 

Comparing expression (|3.29|) with the definition of the decay constant for the nth meson, 

A^ = f n P^n(0), (3.30) 
one can easily conclude that in the chiral limit 



N c 

fn = 8nnV—, (3.31) 



that coincides with relation ( |2.80|) . 
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Figure 8: Graphical representation of the quark-quark scattering amplitude. 



3.4 Quark-quark scattering amplitude 

In this subsection we come to calculation of one of the most fundamental objects in the the- 
ory — the quark-quark scattering amplitude. If known, this object allows one to define the 
dressed current vertices and thus to investigate such properties of the theory as Ward iden- 
tities, current conservation laws, etc. Equation for this object is given in the diagrammatic 
form in Fig.|] and reads 



r aS,-y/3( 7 p 



V (7o), (7o) fem (t% (t% 



(3.32) 



where Greek and Latin letters stand for the colour and the spinor indices, respectively. 
Then one can separate the colour structure of T, 



1 ik,lm — _/y 7 ik > lm > 

and it is also convenient to introduce a new function T defined as 

f 8k,ir(q», k, = S st (q)T tkt i n (q, k, P^S nr (q - P). 



(3.33) 



(3.34) 



Due to the instantaneous type of the interaction induced by the two-dimensional gluon, 
it is also useful to define the amplitude $ as an integral of T over go, 



Then the equation for the new function is 



^ak,ib{p,k,P^ 



(2vr) 2 7 
(p — k) 2 



(7o)ii(7o)fcm 



i J ^-Sai(p)S mb (p - P) 



(3.35) 



(3.36) 



+7 



dq 



:(7o)is(7o)r m $ s fc,/r(<?, k, P^ 



i [ ^-S ai (p)S mb {p - P) 

J Z7T 



(p- q) 2 

The object in the square brackets can be easily integrated out using formula ( |A.1| ) and 
the result reads 

' ' 2^"(P)*«*<P - P)= E(p) + E(P-p)-P 
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+ 



(A_(p)7o)a i (A+(p-P)7o)mb 



(3.37) 



to be 



E(p) + E(P-p) + P 
The simplest way to proceed further is to guess the general structure of the amplitude <£> 



Po ~ PS 



(3.38) 



n=— oo 

so that after some algebraic transformations with the use of the matrix bound-state equation 
Dl) one arrives at 



E xUKP) 



A + (p)$*(p, P)A_(P - p) A_(p)$ n (p, P)A+(P - p) 
E(p) + E(P -p)-P Q E(p) + E(P -p) + P 



(3.39) 



(2tt) 2 7 
(p — k) 2 



(A + (p)) ai (A_(P - p)) kb (A_(p)) aZ (A + (P - p)) 



kb 



E(p) + E(P -p)-P E(p) + E(P -p) + P 



We parametrize the function x n {Pi P) i n t ne form: 

X «(p,P)=27r 7 |dg^ 



(p - q) 2 



(3.40) 



with ip n (p,P) being a new unknown function. Substituting (|3.40|) into (|3.39| ), we arrive at 
two simple equations defining i() n (p,P), 



E S P k n (?> P)n(P, P)} = ±2tt5( P - g). 



(3.41) 



With the help of the formulae from Appendix A the following solution to these equations 
can be found: 

i/j n (p, P) = sign(n)$ n+ G», P), (3.42) 

and, hence, 

X n (p, P) = 27r 7 sign(n) J dq ^^J , (3.43) 
that gives for the quark-quark scattering amplitude f35fl: 



^ik,lmiPi k, P^ 



27TZ7 

(p - k) 



r(7o)«(7< 



O/fcm 



*7 2 E 



Po~ Pfi 



{T n (p,P)) im {T n {k,P)) k i 



+ ^7 2 E 



1 



-( 7o r n (p - P , p)7o) im (7 r„(P - k, p) 7o ; 



(3.44) 



n Po + Po 

where the sum over mesons J2 n counters only positive excitation numbers with n > 0. 

A couple of comments concerning the solution ( |3.44|) is in order. First of all, note that, 
once the wave function $ n contains only symmetric matrices 70 and 75, then = and, 
hence, 

V ik ,i m {p, k, Pfj) = T mlM (p, k, Pfj). (3.45) 

The other comment concerns the form of the solution ( |3.44j) . One can see that the quark- 
quark scattering goes through the formation of a one-meson intermediate state and that the 
last two terms in ( |3.44j ) give nothing but the sum over the full mesonic propagators. 
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Such way, we end with an effective diagrammatic technique which involves the dressed 
quark propagator (|3.1| ) (see ( |A. 1| ) for its ultimate form), the dressed meson-quark-antiquark 
amplitude (r for the incoming and T for the outgoing mesons given by (|3.8| )), the dressed 
quark-antiquark scattering amplitude ( |3.44| ) and the constant —ij/y/Nc which is to be 
inserted into every vertex where a meson couples to a quark-antiquark pair. In addition, 
each quark loop brings two extra factors, the standard fermionic (-1) and Nc from the colour 
trace. One can use these ingredients as bricks for building any hadronic process in the theory. 



3.5 Vector and axial-vector currents conservation. Ward identi- 
ties. 

As stated in the previous subsection, with the quark-quark scattering amplitude in hands 
we are in the position to study properties of currents in the 't Hooft model. Let us prove the 
currents conservation first. For the vector current we are interested in the matrix element 



V* 1 (P) = (n\q 1 ,q\M,P), 



(3.46) 



depicted in Fig.[7|(a). 

On writing the corresponding matrix element and performing the integration over the 
energy p , one finds: 



V? = vyyfih J 7^k S PlS(p ~ Ph,S(p)T M (p, P)] = ^J £^p[7^m(p, P)]. (3.47) 



dp 



(2 

It is easy to check that, multiplying the bound-state equation ( |3.6| ) by •yoy/Nc/Zn, taking 
trace over the spinor indices, and integrating all terms of the resulting equation over the 
momentum p, one arrives at the relation: 



dp 



'N c I ^Sp[ l0 $ M (p,P)} 



dp 



' X r I gSp[ 7 i$Af(p,P)] 



0. 



Combining it with the definition (|3.47 ), one finds that the vector current is conserved, 



P M V M _ py 

Similarly, defining the axial- vector current as 



M 



Af?(P) = (n\q 1 , l5 q\M,P), 



(3.48) 
1, 

(3.49) 
(3.50) 
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Figure 10: Graphical representation of the dressed vector current- quark- ant iquark vertex. 



i.e., using the diagram (b) in Fig.|7] and performing the same steps concerning the bound- 
state equation as before, but with the evident change 70 — > 71 in the multiplier, one can 
derive the axial-vector current divergency in the form: 



pM A M 



PA 



M 



-2m\ Nr 



dk 
2^ 



Sp[7s$ 



M 



(3.51) 



which turns into the axial-vector-current conservation law in the chiral limit. 

It is instructive to see how the same relations appear in the Hamiltonian approach. We 
shall concentrate only on the vector current conservation law as a similar analysis for the 
axial-vector current is straightforward then. 

We start from the definition of the vector current, 



and reformulate it in terms of mesonic creation and annihilation operators,^ 



(3.52) 



Jo(x) = 20Vc [ ^ Px Y,e- lP °" x °m n (P) f f-g n {p, P)/o(p, P) + h.c, (3.53) 

v J III „ J ZTT 



dp 



2?r 
dP 



J{x) = 2^N~ C J ^e^^e- lP o"-o mn(P) J ^-f n (p,P)g (p,P) + h.c. 



(3.54) 



where f n (p,P) and g n (p,P) are introduced in (|2.86| ), whereas f and g are other notations 
for the pion wave functions /„- and g^. 

Then calculating the corresponding matrix element explicitly one finds: 



P M V M _ pyM 



i(n\dMo)\M,p) 



(3.55) 



2JN„ I ^ 



2ti 



P M 9m( P , P)fo(p, P) ~ Pf M (p, P)9o(p, P) 



0. 



where the r.h.s. of this equation vanishes due to the bound-state equation ( |2.40| ) or (|3 . 6|) . 

Now we turn to the investigation of the mesonic form-factors defined by the diagram 
depicted in Fig.[| but we need to know the dressed current-quark-antiquark vertices first. 

Let us start from the vector current. The corresponding expression for the diagrams 
given in Fig.|II] reads 



StoWtMStob - P)T%tt(q,p, P) = (3.56) 



11 Note that we have restored the explicit dependence of operators m n on time in the form m n (xo, P) 
e —iP x °rn n (P) with Pq being the energy of the nth mesonic state moving with the total momentum P. 
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where, as before, the colour indices are denoted by the Greek symbols, whereas Latin ones 
stand for the spinor indices. Substituting the explicit form of the quark-quark scattering 
amplitude T s c kfd (q,P, P) from ( |3.44j ) into ( |3.56| ), one finds 



Fn{p,P) f dq 



Po - PS 



2ty 



Sp[^ n (q,P)} 



7 r n (P -p,P)7o f dq 



J ^ S P^n(P-Q,P)]- 



(3.57) 



Po + P^ 

After tedious but straightforward calculations using the explicit form of the matrix 
bound-state equation ( |3.6|) and the matrix wave function (|3.7|) the following relation can 
be proved: 

- iP^ip, P) = S-\p) - S-'ip - P), (3.58) 

which is nothing but the vector Ward identity for the 't Hooft model (similar expression in 
the light-cone gauge was derived in [|J). 

Now it is a simple task to substitute the solution (|3.58|) into the matrix element written 
for the diagram in Fig.[5] with the vector current instead of the curly line and to arrive at the 
vector current conservation law, 



Q lt (M,P\v lt \M',P') = 0, Q» = P» 



P'. 

A 4 



(3.59) 



Similar calculations give the following results for the axial-vector current in the chiral 
limit |35j , 



T n (p,P) f dq 



Po 

7o r„,(p-p,p)7o 



Sp[l^n(q,P)} 



P n n J 2vr 



^Sp[ ltll ^{P-q,P)l 

Z7T 



and, finally, 



Po + PS 

iP^ip, P) = 5- 1 (p) 75 + l,S-\p - P), 

Pu-P' 



(3.60) 
(3.61) 

(3.62) 



Q,(M,P\ atl \M',P') = 0, ()„ . jr 

Relation ( |3.61|) plays the role of the axial-vector Ward identity. 

In conclusion let us notice that one could arrive at the same results using the Hamiltonian 
approach but at the price of a much more complicated algebra and a much less transparent 
interpretation of the results in terms of Feynman diagrams. 



3.6 The pionic vertex 

The general structure of the mesonic vertex can be considerably simplified in case of the 
pion, since the explicit form of the pionic wave function is known. Indeed, substituting the 
solution ( |2.45| ) into the matrix form ( |3.7| ) and then into the definition ( |3.8| ), one easily finds 



r w (p, p) = s-\p)(\ + 75) - (i - i,)s~\p - p). 



(3.63) 

It is also instructive to compare formula ( |3.63| ) with the Ward identities for the vector 
and axial- vector currents derived in the previous subsection, equations ( |3.58j ) and Q3.61]) , 
respectively. As a result, one finds the following relation for the pionic vertex [[35]: 



r n (p, p) = -iP^ip, p) - iP^{p p) 



(3.64) 
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Note that it is not surprise that the pion couples not only to the axial-vector, but to the 
vector current as well. The reason is that in the two-dimensional theory the axial- vector and 
the vector currents are dual to one another, 

Jl{x)=e tlv J v {x), (3.65) 

where is the totally-antisymmetric Levy-Civita tensor in two dimensions. 



4 Strong decays 

This section is devoted to investigation of hadronic processes in the 't Hooft model at the 
example of the decay A — > B + C. In subsections 4.1 and 4.2 we derive the amplitude of such 
a decay using the Hamiltonian and matrix approach, respectively. We discuss its properties 
and correspondence with nonrelativistic models. Subsection 4.3 is devoted to derivation and 
justification of the two-dimensional Adler selfconsistency condition ("Adler zero"). 



4.1 Suppressed terms in the Hamiltonian 



In section 2 we developed the Hamiltonian approach to QCD2 and diagonalized the Hamil- 
tonian of the model in the mesonic sector. Now let us turn to corrections to the Hamiltonian 
(|2.39|) suppressed by powers of Nc in the denominator. The leading correction is of order 
0{1/ \JNc) and it comes from the terms containing the products MB, MD, M + B, and 
M + D of the operators introduced in (|2.25| ). 



AH 



7 



dpdkdqdQ 6(p) - 9(k) 6(Q - p) - 6(Q - k) 
cos sin — 



N C J (2n) 3 (p - k)' 



x 



M + (p,p- Q)M + {k - Q, q)M{k, q) + M + {p - Q,p)M + (q, k - Q)M(q, k) 
-M + (q,p)M(q,p- Q)M(k -Q,k)- M + (p, q)M{p - Q, q)M{k, k - Q) 



(4.1) 



where M and M + can be defined through the meson creation and annihilation operators 
(see also ( |2.35 )), 



M+{p,k) 



M(p, k) 



E m n( k —P^+ik, k-p) - m n (p - k)<p%{p,p — k) 



n=0 



(4.2) 



E m n {k —p)v+(h, k-p) - m+(p - k)ip\{p,p - k) 

n=0 L 



As easily seen from (|4.1|), this correction describes vertices with three mesons involved. 



4.2 A — ► B + C decay amplitude 

Now we are in the position to study the hadronic processes in the 't Hooft model. Strong 
decays A — > B + C are of most interest for us p5| . As mentioned before, in the leading order 
in Nc the 't Hooft model describes free mesons (see ( |2.39| )), whereas the interaction between 
them is hidden in the suppressed terms partially restored in the previous subsection. Thus 
we expect the amplitude M(A — > B + C) to be of order l/y/N c , whereas M(A + B — > 
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C + D) ~ 1/Nc- In this subsection we study the influence of the backward motion described 
by the tp_ component of the mesonic wave function on the above mentioned amplitude. 

The Hamiltonian approach developed before gives the most straightforward way to calcu- 
late the amplitude of the strong decay since one just needs to evaluate the following matrix 
element: 

M(A^ B + C) = (B(P B )C(P C )\H + AH\A(P A )) = (B(P B )C{P C )\AH\A(P A )). (4.3) 

With the help of the explicit form of the operator AH given in ( |4.1| ) one easily finds the 
general form of the amplitude in terms of mesonic wave functions in the rest frame of the 
decaying particle A {Pa = 0, Pb = —Pc = p) to be (see also f36| where a six-term decay 
amplitude is discussed): 



/ ( -1)2 {-<P-( k +P> V- ( k + P> °) t c (-P' k )v+(^ -P) + S (-P> 9» k )^-c] 



M(A -> B + C)= (4.4) 

7 r dkdq 
m^J (q-kY 

~V+ (k + p, 0)v?+(&, -p) [c(p, q + p,k+ p)ip+(q + p,p) + s(p,q + p,k + p)(p*(q + p, p)} 
-f^(k, -p)<p*(k + p, p) [s(0, q + p,k + p)p+(q + P,0) + c(0, q + p,k + p)(f£(q + p, 0)] 
+<p ( L(k, -p)<f+(k + p, p) [c(0, q + p,k+ p)<p$(q + p,0) + s(0, q + p,k+ p)^(q + p, 0)] 
+ip A (k +p, 0)ip^(k, —p)[s(p, q + p, k + p)<p+(q + p,p) + c(p, q + p,k + p)<p_(q + p,p)\ 
+<fi(k + P, 0)<P+(k + p, p) [s(-p, q, k)ip%(q, ~P) + c(-p, q, k)if^(q, -p)] } 

+ (B <-> C,p <-> -p), 

where 

c(p , q , k) = C0B m^M sm »(P-k)-e(p- q ) ) 
s(f) , k) = am »(*)-»(«) cos Q(P-k)-o ip -, K 

In spite of its frightening appearance, the amplitude ( |4.4j) has a very simple structure. 
Indeed, it contains six terms, i.e., three times more than one could naively expected and 
this is entirely due to the presence of the if- component in the mesonic wave function. 
If one neglects the backward motion contributions in the amplitude (|4.4j ) and inserts the 
nonrelativistic values of the angle 9 (cos8(p) = 1, sm6(p) = p/m), then it reproduces the 
standard quark-model decay amplitude due to the OGE Coulomb interaction |37| adapted 
to the two-dimensional case. It is clear, however, that the substitution of the nonrelativistic 
angle is not justified for kinematical reasons. 

One can arrive at exactly the same expression for the amplitude using the matrix ap- 
proach and writing the matrix element for the diagrams depicted in Fig.p"T|, 

M(A -> B + C) = 

ry 3 f d 2 k 



Sp[T A (k + P B ,P A )S(k-P c )Tc(k,Pc)S(k)T B (k + PB,PB)S(k + P B )] (4.5) 

f N c J {2nr 

+ {B^C). 

In the meantime, the matrix approach proves more convenient in studies of the decay 
amplitude (|4.4j) , (|4.5j ), so we stick with it in the next subsection considering pions in the 
final state. 
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Figure 11: Graphical representation of the amplitude for the decay A — > B + C . 



4.3 Adler selfconsistency condition 

In this subsection we have one more careful look at the pions — namely, at their role in 
the hadronic processes. We remind the reader that these are pions to suffer most drastically 
from the presence of the backward motion of the qq pair described by the ip_ component of 
the wave function. On the other hand, hadronic processes with the pions in the final state 
are much better investigated experimentally, so that any theoretical hint as to how the chiral 
nature of the pion affects hadronic decays is of paramount importance. The 't Hooft model 
for two-dimensional QCD is a source of such hints. 

Thus let the meson B be the pion. We substitute the explicit form of the pionic vertex 
fl3.63| ) into expression ( |4.5| ) and after simple algebraic exercises arrive at a striking conclusion 
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M(A ^n + C) = (4.6) 

in the chiral limit. Note that this result could be anticipated in view of the identification 
Q3.64 ) and the current conservation laws (|3.59| ) and ( |3.62| ). The condition ( |4.6j ) is nothing but 



the two-dimensional analogue of the celebrated Adler selfconsistency condition for amplitudes 



with soft pions involved | 38fl . It is not surprise, that it holds true for any value of the pion 
momentum as, in view of the dimensionlessness of the pion decay constant f n , there is no 
soft scale in the model which could play the role of the edge for this condition. 

Let us conclude with a couple of comments concerning the condition ( }4.6| ). First of all, 
the pion turns out sterile, at least in the sub leading order in Nc- On the other hand, applying 
the above qualitative analysis concerning the dimension of f^, one can extend the condition 



( |4.6|) to the case of any hadronic process with pions involved, so that the pion is completely 
decoupled from the spectrum in the chiral limit in all orders in Nc- 

The formal reason for the condition (|4.6j ) is a totally destructive interference between the 
ip + and the components of the pionic wave function. They were found to be of the same 
order of magnitude and, hence, all six terms of the amplitude (|4.4| ) are equally important 
in establishing the condition ( ^4.6| ). The latter observation automatically invalidates any 
attempts to describe the pion in the framework of a constituent quark model when the y?_ 
component is completely lost. It seems quite unprobable that the above drawback of the 
quark models in QCD 2 , as well as in QCD 4 , could be cured by simple prescriptions, like 



multiplying the naive 2-term amplitude by extra "magic" factors |36[ or whatever 
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5 Conclusions and prospectives 



Phenomenological successes of quark models do tell us that the constituent quarks are the 
correct degrees of freedom in the nonperturbative domain. In these models the lowest 1 Sq qq 
state exists on the same footing as other mesons. As it was already mentioned, there are no 
direct indications that the confinement and the chiral symmetry breaking are interrelated 
in case of QCD 4 . Nevertheless, if such a scenario does not take place, then one easily runs 
into troubles with constituent models: two pions exist, one is a quark bound state, and the 
other one is responsible for the chiral symmetry breaking. Of course, a roundabout way 
is to disregard quark models completely. However, it is more economical to organize the 
confinement and the chiral symmetry breaking due to one and the same mechanism. 

Developing such an approach for QCD 4 is not a straightforward task. Models |16[ 
incorporate the main ingredients, the gap equation similar to ( |2.13|) and the Bethe-Salpeter 
equation similar to ([2.40 ). The existence of the chirally-noninvariant solution of the gap 
equation implies the existence of the Goldstone boson as the lowest qq state. The ip + and 
<f_ components of its wave function are equal to each other in the mesonic rest frame, the 
axial- vector current is conserved in the chiral limit, and all the relations of the current algebra 
are satisfied. In the meantime, the role of the ip_ component of the wave functions for all 



other mesons is less drastic, and parity degeneracy is restored for higher quarkonia [15 



Unfortunately, all these nice features persist at the price of confining interaction em- 
ployed, chosen as the time component of a vector force. First, such a model is not covariant, 
that prevents from proceeding further along the lines described above. In particular, one 
cannot make Lorentz boosts within this model and cannot describe the strong decays. An- 
other drawback is even more important — namely, the interaction is not compatible with 
the area law for the isolated Wilson loop. This point holds true not only for the model [|T5 



where the interaction potential is the three-dimensional oscillator one. It is also so for rather 
sophisticated models with linear confinement too (see, e.g., The reason is that the 

area law yields linear confining potential only for heavy constituents. Besides that, the last, 
but not the least, objection is the lack of the gauge invariance in such kind of models. 



An approach suggested in 17] is rather promising in all these respects. The confin- 



ing interaction employed there is given by a set of gluonic field-strength correlators, which 
produce the area law. These correlators are Lorentz covariant by construction, and gauge 
invariance is preserved too. The latter is ensured by using the generalized Fock-Schwinger 
gauge (Balitsky gauge), which leads to the gluonic correlators explicitly dependent on the 
reference contour and, as a consequence, explicitly translationally noninvariant. A simple 
two-dimensional example of such a correlator is given by the expression ( |2.55| ). The prob- 
lem is very technically involved due to this fact, but at the same time it is very physically 
transparent. Indeed, the interaction of such a kind describes the string which is developed 
between the constituents (for the recent progress in this direction see ^0|). On the other 



hand, the same correlators are responsible for the chiral condensate formation |jl7[ . We 
expect that the quark model followed from such a formalism could be able to describe qq 
bound states including the pion. 

In conclusion, we state once more that a reasonable model in four dimensions is welcome 
in order to find the pionic solution playing the twofold role: being a true qq state it is also 
the Goldstone boson. In our opinion the two-dimensional 't Hooft model gives a brilliant 
example of how this could happen in Nature. 
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6 Appendix A 



In this appendix we collect some formulae useful for the matrix approach. 

The dressed quark Green's function (|3.1|) can be rewritten in a more convenient form if 
the projectors (|2.11|) are introduced, 



A +Obo 



+ 



A_(p)7o 



p - E(p) + ie p + E(p) - is 



(A.l) 



As mentioned before, one should be careful with the sign of the dispersive law E(p) 
and keep the product eE(p), when combining the two fractions in ( |A.1| ) together 0. The 
propagator (|A.1[) contains all radiative corrections and should be assigned to the internal 
quark lines (fat lines in diagrams). 

We find it useful to split the matrix wave function ( |3.7|) into two pieces by means of the 
projectors ( |2.11| ) ,p| 

$„(p, P) = A»$(p, P)A^(P - p), $(p, P) = %(p, P) + § n (p, P), (A.2) 
-, or, vice versa, rj = — , fj = +. Then the completeness condition 



where either 77 = +, 77 
for $'s reads 



£ Sp \$l + (p, P)$ n + (q, P) - ^(p, P)$ n _(q, P)l = 2n5(p - q) 



(A.3) 



which can be easily derived, using the corresponding properties of the functions (p± and the 
relation 



Sp *f(p,P)$%(q,P) =<p%p,P)<f%{q,P) 



9(p) - 9(g) 9(P- p )-6(P-q) 
o m > cos cos 



. 9(p) - 9(g) . 6(P-p)-6(P-g) 
- d V fj> sin sin 



(A.4) 



We also give two formulae useful for the derivation of the quark-quark scattering ampli- 
tude, 

$-"(p,P) = ^ + (P-p,P), (A.5) 



Sp [$"+(P - p, P)$™(P - q, P)j = r, n r, m Sp [$" + (p, P)$™{q, P) 



(A.6) 



2 In this appendix we suppress the index numerating the mesonic states, keeping it only when necessary. 
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where rj and fj are defined above and r\ n is connected to the spatial parity of the state (see 
equation ( |2.44|) ). 

The next two formulae come out directly from the bound-state equation Q2.40| ), 



A„(p)7or(p, P)A^(p - P)7o = - [E(p) + E{P -p)- V P ] A v (p)$(p, P)K n {p - P) (A.7) 

7 

= - [E(p) + E(P -p)- V P ] %(p, P), 

K,(p - PhoT(p, P)HPho = - [E(p) + E(P -p) + rjPo] A v (p - P) 7 o$ P)^)^ 

7 



- [E{p) + E{P -p) + V P ] 7o $+(p, P) 7o , 



(A. 



whereas with the help of relations (|2.43 ) and ( 2.44 ) one finds 



T n (P -p,P) = ^ 7o f _ n (p, P) 7o . (A.9) 
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